Heating in collisions of solids: possible application to impact craters 
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Abstract 

I Due to the importance of coUisions and impacts in early phases of 

04 the evolution of the planetary system, it is interesting to estimate the 

heating of a solid target due to an impact in it . A physically simple 
calculation of the temperature to which a solid target heats up after 
the impact of a projectile with mass m and speed v is performed, and 
possibilities for the application of this result in planetology are pointed 
out0 
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1 INTRODUCTION 



Colliding macroscopic bodies exchange momentum. In inellastic collisions, a 
part of the momentum is transformed into heat. 

The aim of this paper is to analyze to some extent the thermal effects in a 
collision of solid bodies of different masses, and to estimate the temperature 
to which the more massive of the colliding bodies (the target) heats up. 
If the temperature of the target rises sufficiently, melting and ultimately 
vaporisation will occur. The melting temperature of a solid can be estimated 
by the so called Lindemann criterion. The value obtained by the Lindemann 
criterion will be compared with the temperature to which the target heats up 
in a collision. As the temperature to which the target heats up depends on 
the kinetic energy of the impactor this calculation will provide an estimate 
whether or not in an impact the target melts, and as a consequence can it be 
analyzed by applying solid state physics or not. 

If the mass of the target is sufficiently higher than the mass of the im- 
pactor, the impact will produce virtually no change in the momentum of 
the target, while the impactor will almost certainly be destroyed. The ki- 
netic energy of the impactor will be spent on heating of the target and on 
introducing changes into its structure. 

The following section contains the calculation of the temperature to which 
a target in a collision heats up in the case of an inellastic collision. The 
calculation is performed for two cases: for the "ideal" case, in which all 
of the kinetic energy of the impactor is converted into heat, and for a more 
realistic case,in which a part of the kinetic energy is used to form a crack 
in the target. The third part shows how to include the equation of state of 
the material of the target in the calculation in one particular example, and 
the last part of the paper points to possible planetological application of the 
calculation presented in the present paper. 

2 CALCULATIONS 
2.1 The ideal case 

The specific heat of a solid at low temperature is given by ([3]) 
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where A;^ is Boltzmann's constant, h is Planck's constant divided by 27r,T 
the absolute temperature and V the speed of sound waves in the material. 
It can be shown in solid state physics that 

(2) 

where P and p denote,respectively, the pressure and density of the material. 
Inserting eq.(2) into eq.(l),it follows that 

Heating up volume V of the material of the target by one degree requires 
the quantity of energy equal to VCy. The temperature to which this volume 
will heat up when impacted by an object having mass m and speed v is in 
the ideal case given by 

where Tq denotes the intial temperature of the target. Practically, the volume 
V physically represents the volume of a crater resulting from the impact. 
The volume is obviously dependent on the shape of the crater. Assume 
that a crater has the shape of a half of a rotating ellipsoid with distinct 
semi axes,denoted by a,b and c. Physically,a and b are the semi axes of 
the "opening" of the crater, and c denotes its depth. The volume is = 
(2/3)7ra6c,which leads to the final expression for the temperature to which 
this volume of the target heats up as a result of the impact 

Inserting eq.(2) into eq.(5), it finally follows that 

Is this heating sufficient or not to lead to melting and even vaporisation 
of the material of the target? This can be detemined by applying the so 
called Lindemann criterion which states that ([2]) 

Tm = TU^f"" exp(27o(l - (po/p))) (7) 
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In this expression Tm denotes the melting temperature of a material at mass 
density p, T^o is the melting temperature at density po-The symbol 70 de- 
notes the value of the Gruneiscn parameter of the material at density po-The 
Gruneisen parameter is defined as 

oiKt , . 

where a is the thermal expansion coefficient and Kt the isothermal bulk 
modulus of the material. Comparing eqs.(6) and (7) it follows that Ti < 
if 

To 



exp(27o(l - (po/p))) X [1 
x(^)-^/^exp(-27o(l-(po/p))] (9) 

where denotes the kinetic energy of the impactor in the moment of impact. 

In practical terms,the fulfiUement of this condition means that as a con- 
sequence of an impact, the target heats up at the point of impact, but does 
not melt. A further implication is that such an impact can be analyzed by 
solid state physics. The analytical expression of the derivative ^ depends 
on the form of the equation of state of the material of the target. 

2.2 The "real" case 

If the relative velocity of two colliding solid bodies is high enough, and if 
the material strength of the target is small enough, it can be expected that a 
fracture will occur in the target as a consequence of the collision. A standard 
result for the stress necessary for the occurence of a fracture is given by (for 
example, [7]) 

--^(^)"^ (-) 

where E is Young's modulus of the material,x is the surface energy,T the 
radius of curvature of the crack,ro the interatomic spacing at which the stress 
becomes zero and a denotes the length of the crack. A simple analysis shows 
that (7c has the dimensions of pressure, which means that the stress multiplied 
by a volume has the dimensions of energy. 



3 



This means that when a fracture forms in a target as a result of the 
impact the energy ballance has the following form: 



^mv^ - acV = CvV{T, - Tq) (11) 
It follows that in this case the target heats up to a temperature given by 

Ti = To + — x(^-ac) (12) 

Inserting eqs.(2),(3) and the expression for the volume into eq.(12), simple 
algebra could give the result for the temperature Ti to which the target heats 
up. 

A physically more interesting approach to the same problem is given 
by dynamic quantized fracture mechanics {DQFM for short) (for example 
[5]). The basic difference between DQFM and the usual approach used in 
material science is that DQFM introduces geometry in studies of scaling 
laws in material science (for example [1]). Considering that the occurence of 
a fracture in a material is a sign of its failure, it can be shown in DQFM 
that the stress necessary for the occurence of a failure is given by ([5]): 

where Kjc denotes the fracture toughness,po is the radius of the crack of 
length /o and q is the length of the so called fracture quantum. In this case 
the energy balance is determined by eq.(ll), with ct/ instead of ac, and the 
temperature to which the target heats up is given by eq.(12) with the same 
replacement. 

The final result is 

In this way we have obtained an expression for the temperature to which a 
target heats up after the impact in the case in which a fracture forms in it. 
Does it or does it not melt as a result of the impact can again be determined 
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by comparing Ti with the resuh of Lindemann's criterion, expressed by eq.(7). 
It turns out that the condition for " non- melting" is given by 



E, 



x[l-Ax(^)-/3 

J-mO P 



exp(-27o(l - (po/p))] 



(15) 



3 The influence of the equation of state 

Several equations in the preceeding section contain the derivative which 
implies that their application demands the knowledge of the equation of 
state (EOS) of the material of the target. Generally speaking,the EOS is an 
equation of the form P = f{p,T) where all the symbols have their standard 
meanings, and / is some function. The choice of the EOS appropriate for 
any given material is a complex task. 

As an example, in the following the so called Birch-Murnaghan EOS will 
be used. This equation has the form [6] 
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(16) 



where Bq = —V{^)t is the bulk modulus of the material and B'q = (§^)t 
is its pressure derivative. The symbols Vq and V denote the volume of a 
specimen under consideration at the initial value of the pressure Pq and at 
some arbitrary value P. It can be shown that 



and 



3B;)p„p(^)'/» + 5(3B;-16)p?(-!i)2/»l 
Po Po 



^ 1 p, 



0x3/2 



[27{B'^ - 4)p2 + 14(14 - 



Po Po 



(17) 



(18) 



5 



The last expression can be simplified by assuming that B'q = 0. It follows 
that in this case the temperature to which the target gets heated as a result 
of an impact and a fracture forms in it is given by 



The calculations discussed in this paper are mathematically simple, but phys- 
ically they have considerable interest for planetary science. The main re- 
sults, expressing the temperature to which a massive target heats up as a 
result of the impact of a projectile of smaller mass are given by eqs.(6) and 
(16). The calculations have been performed for two distinct cases: when the 
impact has no consequences on the structure of the target, and when it is so 
strong that a fracture occurs in the body of the target. 

There is a distinct difference in the behaviour of the results obtained in 
the two cases. In what we have called "the ideal case", the temperature of the 
target always changes as a result of the impact, as dP/dp is always different 
from zero. No heating effect would occur only in the special case dP/ dp = 
which is physically unrealistic, because a material in which pressure and 
density are not related in some way does not exist. 

An interesting result has been obtained in the "real case". Namely, using 
stadard material science would have given a result which would have taken 
into account only the parameters of the material of the target. Instead of 
such an approach, we have used the DQFM approach, which takes into 
consideration the material parameters, but also the geometry of the problem. 
Note that the result contains the length and radius of the crack, as well as 
the length of the fracture quantum. 

A possible extension of the calculations discussed here could be the inclu- 
sion of the presence of magnetic fields. Recent experiments in material science 
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(19) 



A similar result could be obtained for any other form of the EOS. 



4 Discussion and conclusions 
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have shown that the presence of low magnetic fields can lead to strenghtening 
of materials, [4]. This means that their failure strength is increased. 

In this case the heating effect does not always occur. It is clear that for 
the particular value of the kinetic energy of the impactor, given by 

Ek = afV (20) 

Ti = To in eq. (12), which means that in that case the target does not heat up 
as a consequence of the impact. Such a result may seem strange, but it is 
a consequence of the partition of the kinetic energy of the impactor on the 
heating of the target and formation of a fracture in it. 

Results of this paper have obvious applications in planetary science. It is 
known that impacts and collisions between various bodies have been occuring 
throughout the history of our planetary system. One of the consequences of 
these events are the impact craters on the Earth, Moon, Mars, (at least) some 
asteroids (for example on Eros, see [8]) and satellites. It is expectable that 
the spots on the surfaces of these objects in which impacts occur get heated. 
The question which arises is whether or not such events can be analyzed by 
using sohd state physics. Namely, if the material of the target(s) heats up 
but does not melt or vaporize, an impact can be analyzed by using laws of 
solid state physics. However, if melting occurs, solid state physics is applicable 
only up to the temperature of melting. Details will be discussed elsewhere. 
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